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The sensitivity of the Polyakov Nambu-Jona-Lasinio (PNJL) model as an effective theory of quark 
dynamics to the chiral symmetry has been utilized in studying the QCD phase-diagram. Also, the 
Poyakov linear sigma-model (PLSM), in which information about the confining glue sector of the 
theory was included through the Polyakov-loop potential. The latter is to be extracted from pure 
Yang-Mills lattice simulations. Furthermore, from the quasi-particle model (QPM) as an excellent 
approach to reproduce the lattice QCD calculations, the gluonic sector of QPM is integrated to LSM 
in order to reproduce recent lattice calculations. The hadron resonance gas (HRG) model gives an 
excellent description for the thermal and dense evolution of various thermodynamic quantities and 
characterizes the conditions deriving the chemical freeze-out at finite densities. Therefore, the HRG 
model is used in calculating the thermal and dense dependence of quark-antiquark condensate. 

We review PLSM, QLSM, PNJL and HRG with respect to their descriptions for the chiral phase- 
transition. We analyse the chiral order-parameter M{T), the normalized net-strange condensate 
Aq^s{T) and the chiral phase-diagram and compare the results with the lattice QCD. We conclude 
that PLSM works perfectly in reproducing M{T) and Aq^s{T). The HRG model reproduces Aq^s{T), 
while PNJL and QLSM seem to fail. These features and differences are present in the QCD chiral 
phase-diagram. The PLSM chiral boundary is located in upper band of the lattice QCD calculations 
and agree well with the freeze-out results deduced from the high-energy experiments and the thermal 
models. Also, we find that the chiral temperature calculated from from HRG model is larger than 
that from the PLSM. This is also larger than the freeze-out temperatures calculated in the lattice 
QCD and deduced from experiments and thermal models. The corresponding temperature T and 
chemical potential p sets are very similar to that of PLSM. This might be explained, because the 
chiral T and p are calculated using different order parameters; in HRG vanishing quark-antiquark 
condensate but in PLSM crossing (equalling) chiral condensates and Polyakov loop potentials. The 
latter assumed that the two phase transitions; chiral and deconfinement, take place at the same 
temperature. The earlier deals with the chiral phase-transition independently on the deconfinement 
one. Although the results from the two model PNJL and QLSM keep the same behavior, the chiral 
temperatures deduced from both of them are higher than that of PLSM and HRG. This might be 
interpreted due the very heavy quark masses implemented in both models. 
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I. INTRODUCTION 

At large momentum scale, quantum chromodynamics (QCD) predicts asymptotic freedom [1, 2] or a remark¬ 
able weakening in the running strong coupling. Accordingly, phase transition takes place from hadrons in which 
quarks and gluons are conjectured to remain confined (at low temperature and density) to quark-gluon plasma 
(QGP) [3, 4], in which quarks and gluons become deconfined (at high temperature and density) [5]. Furthermore, 
at low temperature, the QCD chiral symmetry is spontaneously broken; SU{Nf)L x SU{Nf){i —>• SU{Nf)Y. In 
this limit, the chiral condensate remains finite below the critical temperature (Tc). The broken chiral symmetry 
is restored at high temperatures. The finite quark masses explicitly break QCD chiral symmetry. 

Nambu-Jona-Lasinio (NJL) model [6] describes well the hadronic degrees of freedom. Polyakov Nambu-Jona- 
Lasinio (PNJL) model [8-10] takes into consideration the quark dynamics [7] and has been utilized to study 
the QCD phase-diagram [11, 12]. Also, linear-cr model (LSM) [13] can be used in mapping out the QCD 
phase-diagram. 

Many studies have been performed on LSM like 0(4) LSM [13] at vanishing temperature, 0(4) LSM at finite 
temperature [14, 15] and SU{Nf)R x SU{Nf)L LSM for Nf = 2, 3 and 4 quark flavors [16-19]. In order to 
obtain reliable results, extended LSM to PLSM ] can be utilized, in which information about the confining 
glue sector of the theory is included in form of Polyakov-loop potential. The latter can be extracted from pure 
Yang-Mills lattice simulations [20-23]. Also, the Polyakov linear sigma-model (PLSM), and Polyakov quark 
meson model (PQM) [24-26] deliver reliable results. Furthermore, the quasi-particle model (QPM) [27, 28] was 
suggested to reproduce the lattice QCD calculations [29, 30], in which two types of actions are implemented; 
the lattice QCD simulations utilizing the standard Wilson action and the ones with renormalization-improved 
action. 

In the present work, we integrate the gluonic sector of QPM to LSM [31] (QLSM) in order to reproduce the 
recent lattice QCD calculations [32]. In QLSM [31], the Polyakov contributions to the gluonic interactions and 
to the confinement-deconfinement phase transition are entirely excluded. Instead we just add the gluonic part 
of QPM. Therefore, the quark masses should be very heavy. We shall comment on this, later on. In section 
III, we outline the QLSM results. They are similar to that of PNJL. This might be interpreted due the very 
heavy quark masses implemented in both models. Similar approach has been introduced in [33]. The authors 
described in inclusion of gluonic Polyakov loop, which is assumed to generate a large gauge invariance and lead 
to a remarkable modification in hadron thermodynamics. A quite remarkable bridging between PNJL model 
quantum and local Polyakov loop and HRG model has been introduced [34]. A large suppression of the thermal 
effects has been reported and it was concluded that the center symmetry breaking becomes exponentially small 
with increasing the massed of constituent quarks. In other words, the chiral symmetry restoration becomes 
exponentially small with increasing the pion mass. 

The hadron resonance gas (HRG) model gives a good description for the thermal and dense evolution of 
various thermodynamic quantities in the hadronic matter [35-43]. Also, it has been successfully utilized to 
characterizing the conditions deriving the chemical freeze-out at finite densities [44]. In light of this, the HRG 
model can be well used in calculating the thermal and dense dependence of quark-antiquark condensate [45]. The 
HRG grand canonical ensemble includes two important features [38]; the kinetic energies and the summation 
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over all degrees of freedom and energies of the resonances. On other hand, it is known that the formation of 
resonances can only be achieved through strong interactions [46]; resonances (fireballs) are composed of further 
resonances (fireballs), which in turn consist of resonances (fireballs) and so on. In other words, the contributions 
of the hadron resonances to the partition function are the same as that of free (non-interacting) particles with 
an effective mass. At temperatures comparable to the resonance half-width, the effective mass approaches 
the physical one [38]. Thus, at high temperatures, the strong interactions are conjectured to be taken into 
consideration through the inclusion of heavy resonances. It is found that the hadron resonances with masses up 
to 2 GeV are representing suitable constituents for the partition function [35-43]. Such a way, the singularity 
expected at the Hagedorn temperature [35, 36] can be avoided and the strong interactions are assumed to be 
taken into consideration. Nevertheless, validity of the HRG model is limited to the temperatures below the 
critical one, T^. 

In the present paper, we review PLSM, QLSM, PNJL and HRG with respect to their descriptions for the 
chiral phase-transition. We analyse the chiral order-parameter M(T), the normalized net-strange condensate 
Ag^s{T) and the chiral phase-diagram and compare the results with the lattice QGD [47-49]. The present work is 
organized as follows. In section II, we introduce the different approaches SU(3) PLSM [50] (section II A), QLSM 
(section IIB), PNJL (section IIG) and HRG (section HD). The corresponding mean field approximations are 
also outlined. Section HI is devoted to the results. The conclusions and outlook shall be given in section IV. 


II. SU(3) EFFECTIVE MODELS 
A. Polyakov Linear Sigma-Model (PLSM) 

As discussed in Ref. [31, 50]. the Lagrangian of LSM with Nf = 3 quark flavors and Nc = 3 (for quarks, 
only) color degrees and with quarks coupled to Polyakov loop dynamics was introduced in [26, 51], 

C^C,mrai-U{(t>,ct>*,T). (1) 

where the chiral part of the Lagrangian of the SU(3);X SU(3)fl symmetric linear sigma-model Lagrangian with 
Nf = 3 is [52, 53] Cchirai = Cq -\- Cm- The first term is the fermionic part, Eq. (2) with a flavor-blind Yukawa 
coupling g of the quarks. The second term is the mesonic contribution, Eq. (3) 

A = gTa{aa+i"t5Tra))tpf, (2) 

/ 

Cm = Tr(a^ 4 >'fa^$-m^ 4 >'f$)-Ai[Tr($^$)] 2 -A 2 Tr($'f$)^-kc[Det( 4 >)-fDet($ 1 ')]-fTr[i 7 ($-k$^)]. (3) 

The summation runs over the three flavors (f=l, 2, 3 for the three quarks u, d, s). The flavor-blind Yukawa 
coupling g should couple the quarks to the mesons. The coupling of the quarks to the Euclidean gauge field 
is given via the covariant derivative [20, 54]. In Eq. (3), $ is a complex 3x3 

matrix which depends on the era and tTq [53], where 7 ^ are Dirac 7 matrices, a a are the scalar mesons and tTq 
are the pseudoscalar mesons. 

$ = Ta(t>a = Ta{cra + IT^a), (4) 

where Ta = Aa/2 with a = 0, • • • ,8 are the nine generators of the U(3) symmetry group and Aa are the eight 
Gell-Mann matrices [13]. The chiral symmetry is explicitly broken through 

H = Taha- (5) 

which is a 3 X 3 matrix with nine parameters ha- Three finite condensates do, do and cfg are possible, because 
the finite values of vacuum expectation of $ and $ are conjectured to carry the vacuum quantum numbers and 
the diagonal components of the explicit symmetry breaking term, ha, where ho 0, ho = 0 and hg 7^ 0 , and 
squared tree level mass of the mesonic fields mf, two possible coupling constants Ai and A 2 , Yukawa coupling 
g and a cubic coupling constant c can be estimated as follows, c = 4807.84 MeV, hi = (120.73)^ MeV^, 
hs = (336.41)3 ^2 ^ _(306.26)2MeV2, Ai = 13.48 and A3 = 46.48 and g = 6.5. 

In order to get a good analysis it is more convenient to convert the condensates cro a-nd erg into a pure 
non-strange ax and strange ay condensates [55] 

f crx\ ^ J_ I' 1 \ 

\<^v j r/3 V 1 y Ws / 


(6) 


4 


The second term in Eq. (1) ,T) represents Polyakov-loop effective potential [20], which agrees well 

with the non-perturbative lattice QCD simulations and should have Z(3) center symmetry as pure gauge QCD 
Lagrangian does [9, 24]. In the present work, we use the potential [/(</), ^*,T) as a polynomial expansion in (j) 
and <j)* [9, 10, 56, 57] 


where 




^2(T) 2 ^3 ^^3 , 

+< 




i>2{T) — ttg + ttl 



+ 02 




(7) 

( 8 ) 


In order to reproduce pure gauge lattice QCD thermodynamics and the behavior of the Polyakov loop as a 
function of temperature, we use the parameters oq = 6.75, oi = —1.95, 02 = 2.625, 03 = —7.44, 63 = 0.75 and 
64 = 7.5. For a much better agreement with the lattice QCD results, the deconhnement temperature Tq in pure 
gauge sector is fixed at 270 MeV. 


1. Polyakov Linear Sigma-Model (PLSM) in Mean Field Approximation 

In thermal equilibrium the grand partition function can be defined by using a path integral over the quark, 
anti-quark and meson helds 


Z 


Trexp[-(7i 


f—u,d,s a 


Vip'D'tpexp 


f{C+ E M/V’/7V/) 

J X jr , . j _ 


f—u,d,s 


(9) 


where f^ = i dt Jyd^x and V is the volume of the system, /i/ is the chemical potential for / = {u,d,s). 
We consider symmetric quark matter and define a uniform blind chemical potential fXf = Pu,d = Ps- Then, we 
evaluate the partition function in the mean field approximation [53, 58]. We can use standard methods [59] in 
order to calculate the integration. This gives the effective potential for the mesons. 

We define the thermodynamic potential density of PLSM as 

= U{<7^,ay) ,T) A- (10) 

Assuming degenerate light quarks, i.e. q = u,d, the quarks and anti-quarks contribution potential is given as 
[51] 


= -2TNq j ^0^{ln[l + 3(</. + cj )*x 

-|-ln[l -I- 3((()* -b X -b 

-2TNs f -^{ln[l-b3((/>-b((-*e-(®»-'')/'^) X 
J [^Trp 

-bln[l -b 3{4>* -b X -b 6“^^®“+'")/^]}, (11) 


where Nq = 2, = 1, and the valence quark and antiquark energy for light and strange quark, Eq = rriq 

and Eg = respectively. Also, as per [55] the light quark sector, Eq. (13), decouples from the strange 

quark sector (m^) and light quark mass rriq gets simplified in this new basis to 


mg 

rus 


2 ’ 
(Jy 

'71‘ 


( 12 ) 

(13) 











5 


The purely mesonic potential is given as 


U{(Jx,cFy) — 2 (a^ + ay) hxUx 2-\/2 


+ ^ 2 ) 0 -^ + ;^(Ai + A2)cry. (14) 


We notice that the sum in Eqs. (7), (11) and (14) give the thermodynamic potential density similar to Eq. 
(10), which has seven parameters m?, hx, hy, Ai, A 2 , c and g, two unknown condensates Ox and Gy and two order 
parameters for the deconfinement (j) and (j)*. The six parameters m^, hx^hy^Xi, A 2 and c are fixed in the vacuum 
by six experimentally known quantities [53]. In order to evaluate the unknown parameters Ox^ Gy, (j) and (jf, 
we minimize the thermodynamic potential, Eq. (10), with respective to Gx, Gy, (f) and (jf, respectively. Doing 
this, we obtain a set of four equations of motion 


dn _ dn _dVt _ dn 

dGx dGy d(j) d(j>* 


^= 0 , 


where min means Gx = Gx, Gy = Gy, (j) = (j) and (jj* = (p* are the global minimum. 


(15) 


m 


XlpLSM — ^ (0; 0; T) T (^a: T ^xGx ^y^y 


2y/2 


^12 2 

Y ^ y 


1 


- 2TN„ 


[ f f {ln[l + 3{cj) + (j)*e X e + e ^^/^] 

J (27r)'^ 


^ ' (2^)3 

+ln[l + 3(^* + X + e“3(W+M)/rj| 

- 2TN. 


(Pp 


(27r)^ 


{ln[l + 3((() + X + e-3(^^»-/^)/^] 


+ln[l + Z{(j)* + ())e“(-®'’+'^)/^) X + e“3(-Eo+A‘)/T’]},^ 

Accordingly, the chiral order parameter can be deduced as 

{tp'lp)pLSM IPs d^lpLSM 


MpLSM = PT-s 


2^4 


dmi 


1 


4" o (2Ai + X2)Gx, + 7(Ai + X2)g, 


(16) 


(17) 


B. Linear Sigma-Model and Qnasi-Particle Sector (QLSM) 

When the Polyakov contributions to the gluonic interactions and to the confinement-deconfinement phase 
transition are entirely excluded, the Lagrangian of LSM with TV/ = 3 quark flavors and W = 3 (for quarks, 
only) color degrees of freedom, where the quarks couple to the Polyakov loop dynamics, has been introduced in 
Ref. [26, 51], 

C^C,Mrai-U{4>,ct>*,T). (18) 

The main original proposal of the present work is the modihcation of Eq. (18) 

C- — Cchiral — p), (19) 

where the chiral part of the Lagrangian Cchirai = Cq + Cm is of SU(3 )lX SU(3)/i symmetry [52, 53]. Instead 
of U{(j), (jf ,T), the gluonic potential p), which is similar to the gluonic sector of the quasi-particle model, 

is inserted, (review Eq. (35)). The Lagrangian with Nf = 3 consists of two parts; fermionic and mesonic 
contributions, Eqs. (2) and (3), respectively. 

Some details about the quasi-particle model are in order. The model gives a good phenomenological descrip¬ 
tion for lattice QCD simulation and treats the interacting massless quarks and gluons as non-interacting massive 
quasi-particles [60]. The corresponding degrees of freedom are treated in a similar way as the electrons in con¬ 
densed matter theory [61], i.e. the interaction with the medium provides the quasi-particles with dynamical 

masses. Consequently, most of the interactions can be taken into account. When confronting it to the lattice 

QCD calculations, the free parameters can be fixed. The pressure at finite T and p is given as 

P= '^Pi-B{T,p), = dk [f+(k) + f-(k)] , (20) 
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where the function B (T, stands for bag pressure at finite T and ^ which can be determined by thermodynamic 
self-consistency and dp/dIVa = 0 ; the stability of p with respect to the self-energies (Ila) and the distribution 
function for bosons and fermions, ± respectively is given as 

ft[k) = — 

exp 


Ei{k)TtJ. \ 


± 1 


The quasi-particle dispersion relation can be approximated by the asymptotic mass shell expression near the 
light cone [27, 28], 

Ef{k) = +m^,{T,p) = e + Il,{k;T,p) + {x,Tf, ( 22 ) 

where Ili(k;T, p) is the self-energy at finite T and p and xf is a factor taking into account the mass scaling 
as used in the lattice QCD simulations. In other words, xf was useful when the lattice QCD simulations have 
been performed with quark masses heaver than the physical ones. In the present work, the gluon self-energies 
Ilg{k-,T,p) are relevant [62] 


Nf 


- ■ f 


2 3 V- 2 \ 


Ug{k]T,p) = 

where the effective coupling G at vanishing chemical potential is given as, 

C? 2 (T) = 


Gloo^in 


T>Tc, 


(23) 


Gloop(T) + b{i-^), r<r,. 
And the two-loop effective coupling G 2 ioop(^) reads [27] 


(24) 


^2 (T) - [l _ oZZZZl 

^ /3oln(e) [ /32 ln(e) J 


e = A 


T-T, 

Tr 


(25) 


and Tg is a regulator at T^. For The parameter A is used to adjust the scale as found in lattice QCD simulations. 
These two parameters are not very crucial in the present calculations. The regulator and scale are controlled by 
the condensates and ay) and the order parameters (</) and </)*), which are given as function of temperature 
and baryon chemical potential. The fi function [63] depends on the QCD coupling G, P = dG/{d In(A^)), with 
Afj, is the energy scale. It is obvious that the QCD coupling in Eqs. (24) and (25) is given in T-dependence, 
only. In calculating j3 = dG/{d In(A^)) at finite p it is apparently needed to extend G to be /r-dependent, as 
well. The two-loop perturbation estimation for /3 functions gives 


/3o = -^ (Uric (26) 

/3i = i (34n2-13n/nc + 3^ ) . (27) 

6 V ricJ 


1. Linear Sigma-Model and Quasi-Particle Sector (QLSM) in Mean Field Approximation 


As in section IIA I and Eq. (9), we derive the thermodynamic potential density in the mean field approxi¬ 
mation. This consists of three parts: mesonic and quasi-gluonic potentials in additional to the quark potential 

P) = = u{ag,, ay) + Ug{T, p) + (28) 

• First, the quark potential part [53] 

oo 

/ d^k 

—-^In [I — nqj{T, pf)] -|- In [1 — riqjiT, pf)] 

J=U,a,S g ' 


{ 29 ) 
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It is obvious that is equivalent to 17,,. The occupation quark/aniquarks numbers read, 




( 30 ) 


and antiquarks riqjiT, fif) = nqj{T, respectively. The number of internal quark degrees of freedom 

is denoted by dg = 2 and Nc = 6 (for quarks and antiquarks). The energies are given as 


EqJ — \l k'^ + 




(31) 


with the quark masses m/ which is related to niq and nis for m-, and d- and s-quark, respectively. As 
given, the latter are proportional to the u-fields 


TTlq = g 


rus = g 




(32) 

(33) 


where the Yukawa coupling g = 8.3. The symbols for the chiral condensates, ax and ay for light- and 
strange-quarks, respectively, are kept as in literature. 


Second, the purely mesonic potential part reads 


U{ax,ay) — (tTj, -|- CTj^) hxax hyay ^axay-\- „ 


Ai 


(2Ai-I-A 2 ) cr^ (Ai-I-A 2 )(t: 


2 \~x • -y/ --y- y ^ ' 2 

• Third, the quasi-gluonic potential part is constructed from Eqs. (22), (21) and (20) 


Va = -- 


6 TT^ 


k'^ dk- 


1 


1 


1 


exp — 1 exp 


- 1 


■ (34) 


(35) 


In Eq. (35), the degeneracy factor dg = 8 and two parameters A and Tg, which were given in Eq. (25), 
should be fixed in order to reproduce the lattice QCD calculations. Here, we find that A = 2.0 and 
Tg = 0.0 MeV give excellent results. 

When adding the three potentials given in Eqs. (35), (34) and (29), the thermodynamics and chiral phase 
translation can be analysed. The resulting potential ^Iqlsm can be used to determine the normalized 
net-strange condensate and chiral order parameter, Eq. (55). 


^QLSM = + ay) - hxax - hyay - + ^2)o'x + ^(^1 + ^ 2 ) 0 '.^ 


y 


d r°° 

I dk [21n(l - f-{T,g)) + 21n(l - /+(T,/r)) + ln(l - friT,g)) + ln(l - /+(T,/r))] 


— d, 


J k'^ dk 


(36) 


where, 




ftiT,g) = 




I 


yj +k‘^ ±y, 

e - 'F --h 1 

1 

^9'^ gy {I', +k'^±y 

e T -I-1 


8 ^ 2/9 ^ I oW2^ (l 3(34iV,^-13Ar.iV,+3^)ln(ln^(t^)) ^ 

\'27^h+y2 +OJI Jli (llAr,,- 2 Arj )2 1 n^({ 2 ) 

(IliVc-2 7V/)ln2 (^2) 


1/2 


(37) 

(38) 

(39) 



























^ is function of T and the quark masses should be very heavy. The QLSM results are similar to that 
of PNJL, section III. This might be interpreted due the very heavy quark masses implemented in both 
models. 


C. Polyakov Nambu-Jona-Lassinio (PNJL) Model 


The Lagrangian of PNJL reads [64, 65] 

E 

f=u,d,s f f 

+ Y ^ - gD[dettpfPLtpf'+ det'ijjfPR'iljff] 

a—0,...,S 

-U{(l>[Al^[A],T) (40) 


where the matrices Pl.r = (l± 75)/2 are chiral projectors, l{{cj)[A\, ^[A],T) is the Polyakov loop potential 
(Landau-Ginzburg type potential) and stand for gauge field interactions. The mass of 

a particular flavor is denoted by m/, where / = u,d,s. The two coupling constants go and gs, A“ are 
Gell-Mann matrices [13] and 7 ^ are Dirac 7 matrices. The model is not renormalizable so that we have 
to use three-momentum cut-off regulator A in order to keep quark loops finite. 

The Polyakov loop potential is given by [ 66 ], 


with 


r4 = — 


6 


■4>^) + 


(41) 




hiT) = ao + ai 



-t- 02 




63 and 64 being constants and we choose the following fitting values for the potential parameters, oo = 6.75, 
oi = —1.95, 02 = 2.625, 03 = —7.44, 63 = 0.75 64 = 7.5 and Tq = 187 MeV. These are adjusted to the 
pure gauge lattice data such that the equation of state and the Polyakov-loop expectation values are 
reproduced. 


1. Polyakov Nambu-Jona-Lassinio (PNJL) Model in Mean Field Approximation 


The thermodynamic potential density of PNJL is defined as 

9d 


fl =ll[(j),(j),T]-\-2gs ^ crj --^au(Td(7s - 6 

f=u,d,s f 

(Pp 


tfp 

2— ^ 




Yjj 

/ -^0 


0 (27r)" 

d^p 


0 ( 27 r) 


■ In 


In 


1 -b 3 U) -b (j)e - 'T' - ) e-5^- 


-(Epf+P)\ -(Epf+P) 

1 + 3 ( <)> + ((e -^- - 


EpfOiA - \p\) 

-3(Epy-p)' 

be-- 

- 3 (Epf+p)~ 


(42) 


where Epf = ^p^ -\- Mj is the single quasi-particle energy, = cr^ -b cr^ -b and from isospin symmetry, 

Oq = cr« = (Jd- In the above integrals, the vacuum integral has a cutoff A whereas the medium dependent 
integrals have been extended to infinity. By the self-consistent gap equation, the quark mass can be 
estimated. 
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rriu [MeV] 

ms [MeV] 

A [MeV] 


goJP 

5.5 

134.758 

631.357 

3.664 

74.636 


Tab. I: Parameters of the SU(3) PNJL model. 


where erf = (ipfipf) denotes the chiral condensate of quark with flavor / and other parameters are listed 
out in Tab. I [64, 65]. For isospin symmetry, we define the light and strange-quark masses as 


Ms 

0 , 9d 2 

= ms- 2gs as + — aq, 

(44) 

Mq 

- 0 1 5-D 

— TTig 2ggaq + ^ aq as- 

(45) 


Here, we notice the strong dependence on the cr-fields. 

Now, we have all the PNJL Model parameters except aq,(7s,4> ^-nd </>, which can be estimated from 
minimizing the thermodynamic potential, Eq. (47), with respective to cr^, cts, cj) and (j), respectively. Doing 
this, we obtain a set of four equations of motion 


dVt_dVt_dVl_dVl 

f~v f~v 1 ‘~\~f 1771271 

(jej) (jcj) 


Then, the potential of the PNJL model reads 


^PNJL =l^['P,4>jT] + 2gs gf - 

f — U,d,S 

(Pp 


/ - 

-2Er 

/ “ 


0 

r <pp 

0 ( 2 ^)" 


In 


In 


1 + 3 

1 + 3 


^auCTdCTs - 6E / - 1^) 

^ j: Jo [ZTT) 


- 




-3(Epf-p) 

(j) + (pe 

T 

1 e T 

+ e 

T 

- 

-(Epj.+p) ^ 

\ -(®p/+e) 



(j) + (pe 

T 

1 e T 

+ e 

T 


(46) 


(47) 


Having completed the introduction of both PLSM and PNJL, it is in order now to discuss central Z{3) 
symmetry related to the Polyakov loop. It has been shown that the SU(3) color-singlet has Z{3) symmetry 
through the normalized character in the fundamental representation of SU(3), $( 0 i, 02 ). This becomes 
equivalent to an ensemble of Polyakov loop [67]. Furthermore, it was concluded that 4>(0i, 6 * 2 ) can taken as 
an order parameter for color-confinement to color-deconfinement phase transition, i.e. the center symmetry 
is spontaneously broken at high temperatures. 

Ref [ 68 ] introduced an attempt to resolve some incongruities within NJL and PNJL. It was argued that by 
integrating corresponding extremum conditions, the thermodynamic potential is directly obtained, where 
the integration constant can be fixed from Stefan-Boltzmann law. Keeping the regulator finite at finite 
temperature and chemical potential is the main advantage of this approach. 


D. Hadron Resonance Gas (HRG) Model 


Treating hadron resonances as a free (non-interacting) gas [35-39] is conjectured to give an accurate 
estimation for the thermodynamic pressure below T^. It has been shown that thermodynamics of strongly 
interacting system can also be approximated as an ideal gas composed of hadron resonances with masses 
< 2 GeV [38, 69], i.e. confined QCD matter (hadrons) is well modelled as a non-interacting gas of 
resonances. The grand canonical partition function reads 

^ N-H 

exp T 


ZiT,p,V) = TV 


(48) 
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where T (fj.) is temperature (chemical potential). The Hamiltonian (H) is given as the kinetic energies of 
the relativistic Fermi and Bose particles. 

The main motivation of using H is that 

— it contains all relevant degrees of freedom of confined, strongly interacting QCD matter, 

— it implicitly includes interactions, especially the ones leading to formation of resonances and 

— it gives a quite satisfactory description of the particle production in heavy-ion collisions. 

With these assumptions, the thermodynamics is resulted from single-particle partition functions 


lnZ(T,/r„H) 



i 



1 ± exp 


fli 

T 


dk^ 


(49) 


where £i{k) = (fc^ -|- mf )^/^ is the i—th particle dispersion relation, gi is spin-isospin degeneracy factor 
and ± stands for bosons and fermions, respectively. 

For hadron resonances which not yet measured, experimentally, a parametrization for a total spectral 
weight has been proposed [70] as a recent estimation for Hagedorn mass spectrum [71]. In the present 
work, we merely include known (measured) hadron resonances with mass < 2 GeV. This mass cut-off is 
assumed to define the validity of HRG in modelling the hadronic phase. Resonances with heavier masses 
diverge all thermodynamic quantities at the Hagedorn temperature [35, 36]. 

Very recently, it has been shown that indeed the viral expansion is a reliable way to include hadron 
resonances, because the phase shift is a directly-accessible quantity in experiments [72]. For instance, 
for accurate isospin-averaged observables, the scalar-isoscalar /o( 500 ) (cr meson) resonance and scalar 
Rr*( 800 ) should be not be included in the HRG model. 

The HRG model has been used in calculating the higher-order moments of the particle multiplicity, in 
which a grand canonical partition function of an ideal gas with experimentally-observed states up to a 
certain mass cut-off is utilized [44]. The HRG model has been successfully utilized in characterizing two 
different conditions generating the chemical freeze-out at finite densities, namely constant normalized- 
entropy density s/T^ = 7 [73-76], constant product of kurtosis and variance = 0 [77] and constant 
trace-anomaly (e — 3p)/T'^ = 7/2 [78]. As introduced in [76], the third freeze-out conditions, which is 
characterized by constant s/T^ is accompanied by constant s/n. 

Our HRG model was used to study the possible differences between the behavior of light {qq) = (uu) = (dd) 
and strange (ss) quark-antiquark condensates in hadron phase. The contribution to the pressure due to 
a particle of mass m/j, baryon charge B, isospin J 3 , strangeness S, and degeneracy g is given by 


Ap = 


9 ml 

27r^ 


00 

E 

n—l 


i-v) 


n+1 


exp ( n 


Bp, B — hl^i — Sp,s 
T 


K, 


( rnh '\ 
W T J 


(50) 


where Kn{x) is the modified Bessel function. In hadrons, the isospin is an almost exact symmetry. 

The quark-antiquark condensates are given by the derivative of Eq. (50) with respect to the constituent 
quark masses 


im) 


(ss) 


(99)0 + 

h 


dmh dAp 

dmq dmh ’ 


(■ss)o +E 

h 


dmh dAp 
dms dmh ’ 


( 51 ) 


where {qq)o and (ss)o are light and strange quark-antiquark condensates in vacuum, respectively [45]. It 
was found that at small chemical potential the strange quark-antiquark condensate is larger that the light 
one. At large chemical potential, such difference gradually diminishes. 

Some authors still prefer to take into account repulsive (electromagnetic) van der Waals interactions in 
order to compensate the strong interactions in hadron matter [78]. Accordingly, each resonance constituent 
is allowed to have an eigen-volume. Thus, such total volume should be subtracted from the fireball volume 
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or that of the heat bath. Also, considerable modifications in thermodynamics of hadron gas including 
energy, entropy and number densities are likely. The hard-core radius of hadron nuclei can be related to 
the multiplicity fluctuations. 

About ten years ago, Tawfik derived 5'-matrix for the HRG model [38], which describes the scattering 
processes in the thermodynamical system [79]. Accordingly, Eq. (49) can be written as an expansion of 
the fugacity term 

OO 

lnZ(“*)(E,T,Ai) = lnZ(*'')(E,r,^) + ^a,(r)exp(Ai,/T). (52) 

j /=2 


where ai,{T) are the virial coefficients and the subscript v refers to the order of multiple-particle interac¬ 
tions. 


/ OO o 

dw exp{-e^(u;)/r} y^(2/ -b l)—8i{w). 


(53) 


The sum runs over the spatial waves. The phase shift 5i{w) of two-body inelastic interactions, for instance, 
depends on the resonance half-width T^, spin and mass of produced resonances. 


In Z(“‘)(E,T,/r) 


In n) + 


9r n Tr exp{{-er-{w) + fJ,r)/T} 

/ “It;-- 


(54) 


In Eq. (54), by replacing fj, by —/r, the anti-particles are taken into consideration. For a narrow width 
and/or at low temperature, the virial term decreases so that the non-relativistic ideal partition function 
of hadron resonances with effective masses is obtained. This means that, the resonance contributions 
to the partition function are the same as that of massive free resonances. At temperatures comparable 
to r^, the effective mass approaches the physical one. Thus, we conclude that at high temperatures, the 
strong interactions are taken into consideration via heavy resonances, Eq. (49), i.e. Hagedorn picture. We 
therefore utilise the grand canonical partition function, Eq. (52), without any corrections. 

In order to verify this picture, Tawfik checked the ability of HRG with finite-volumed constituents in 
reproducing lattice QGD thermodynamics [78]. At radius r > 0.2 fm, the disagreement becomes obvious 
and increases with increasing r. At high temperatures, the resulting thermodynamics becomes non¬ 
physical. It was concluded that the excluded volume seems to be practically irrelevant. It has e negligible 
effect, at r < 0.2 fm. On the other hand, a remarkable deviation from the lattice QGD calculations 
appears, especially when the radius r become large. 

In the present work, the chiral parameters, M{T) and Aq_s(T), see section III, are extracted from HRG 
assuming fully and partially chemical non-equilibrium [80]. There is no difference, when yg = I.O and 
when it is allowed to have values different than unity, where 7 g and yg refer to non-equilibrium treatment 
or occupation factors for light and strange quarks, respectively. These two parameters enter Eq. (49) 
after raising them to exponents reflecting the light and strange quarks contents of *-th hadron. They are 
identical to the fugacity factor and therefore are multiplied to the exponential function. 


III. RESULTS 


A systematic comparison between PLSM, PNJL, QLSM and HRG is presented. It intends to calculate 
two chiral quantities, the order parameter M{T) and the normalized net strange non-strange condensate 
Ag^s{T). The results shall be confronted to the lattice QGD simulations [47-49]. The comparison with the 
lattice should signal which model is close to the lattice and on other hand offers differentiation between 
the SU(3) effective models, themselves. 


A. Chiral order-parameter M{T) 


The chiral order-parameter M{T) was originated in lattice techniques [81]. The latter calculates dimen¬ 
sionless quantities in units of lattice spacing rather than physical units. The lattice spacing can then be 
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Fig. 1 : The thermal behaviour of the dimensionless chiral order-parameter, M, calculated as function of temperature 
from the four SU(3) effective models, PLSM (solid curve), QLSM (dotted curve), PNJL (double-dotted curve) and HRG 
(dash-double-dotted curve) and compared with the lattice QCD calculations (solid circles) [47] and (open circles) [48] at 
mi/nis = 0.037. 


converted into the physical units. M(T) relates the light quark condensate to the strange quark mass 


M(T) 


rpi 


(55) 


where iris is the strange quark physical mass which fixed here to 138 MeV in order to get the ratio of 
light and strange quark irrasses mi/ms = 0.037. Also, we notice that the diirrensionless M{T) depends on 
the thermal behavior of the light quark condensate {'tptpi)- Iii lattice QCD, the chiral condensate remains 
finite. But it contains contributions which would diverge in the continuum limit. Therefore, it requires 
renormalization, in particular an additive and multiplicative renormalization. In order to remove - at 
least - the multiplicative renomalization factor, we take into consideration Eq. (55) as a definition for 
the order parameter. The light quark condensate itself can be calculated from the potential, PLSM: Eq. 
(16), QLSM: Eq. (36), PNJL: Eq. (47) and HRG: Eq. (51). Accordingly, we estimate Eq. (55) from the 
four models and then compare them with the lattice QCD calculations. Fig. 1. We find that this chiral 
order-parameter in the SU(3) effective models and first-principle lattice QCD simulations [47, 48] rapidly 
decreases with increasing T. 

Comparing with the lattice QCD [47], the best agreement is found with PLSM, but PLSM underestimates 
the recent lattice QCD [48]. In fact, the lattice calculations [48] lay on top of all curved from the SU(3) 
effective models. This might be originated to the specific configurations of the lattices and the actions 
implemented in each simulations, section III A. The other effective models lay below the two sets of lattice 
calculations. 

The four models PLSM, QLSM, PNJL, HRG and different sets of lattice calculations have different critical 
temperatures. In Tab. H, we list out the critical temperatures corresponding to each order parameter. In 
determining the pseudocritical temperatures, different criteria are implemented. They are not only 
quite unorthodox but they are distinguishable from each other. Further details shall be elaborated in 
Section HI C. 
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Tx [MeV] 

Order Parameter 

PLSM 

164 

crossing of {ox, Oy) and {(f, cj)) 

QLSM 

200 

{cJx, Oy) and largest fluctuation in 1712 / 

PNJL 

217 

crossing of [ox, Oy) and (</>, cj)) 

HRG 

184 

vanishing (i/)'(/))-condensate 

LQCD [47] 

156 at nis/mi = 0.037 

sudden drop in M{T) and Ai^s{T) 

LQCD [48] 

165 - 170 

sudden drop in M{T) and Ai^sjT) 

LQCD [49] 

165 

sudden drop in M{T) and Ai^s{T) 


Tab. II: The pseudocritical temperatures as calculated from PLSM, QLSM, PNJL, HRG and the different sets of 
lattice QCD calculations. 


1. A short comparison between the two sets of lattice QCD calculations 


Refs. [47, 48] presented results for 2 + 1 quark flavors, where all systematics are controlled, the quark 
masses are set to their physical values and the continuum extrapolation is carried out. Larger lattices 
and a Symanzik improved gauge besides a stout-link improved staggered fermion action are implemented. 
Depending on the exact definition of the observables, the remnant of the chiral transition is obtained at 
Tc = 150 MeV. Extending these results, the transition temperature was also determined for small non¬ 
vanishing baryonic chemical potentials. At high temperatures, the lattice pressure is found ~ 30% lower 
than the Stefan-Boltzmann limit. 

Ref. [49] used 2-1-1 quark flavors with physical strange quark mass and almost physical light quark 
masses. The calculations have been performed with two different improved staggered fermion actions, 
the asqtad and p4 actions. Overall, a good agreement between results obtained with these two 0{a^) 
improved staggered fermion discretization schemes is found. At high temperatures, the lattice pressure is 
^ 14% lower than the Stefan-Boltzmann limit. 

From this short comparison, we find that: 

— [47, 48] implement Symanzik improved gauge and stout-link improved staggered fermion action. The 
resulting pressure is found ~ 30% lower than the Stefan-Boltzmann limit. 

— [49] uses improved staggered fermion actions; the asqtad and p4 actions. The resulting pressure is 
^ 14% lower than the Stefan-Boltzmann limit. 


2. Couplings in PLSM 


In the effective models, the parameters, especially the couplings, are very crucial for the outcome of the 
calculations. One of the motivations for the present work is the failure of PLSM [50] in reproducing the 
lattice QCD results [47-49] even with large coupling g. In Ref. [50], g ranges between 6.5 and 10.5. The 
first value was enough to reproduce the lattice QCD calculations, PRD80, 014504 (2009) and PLB730, 
99 (2014). Increasing g to 10.5 does not enable PLSM to reproduce the other lattice simulations [47-49]. 
Furthermore, through fitting with lattice QCD calculations and experiments, the parameters of PLSM 
can be estimated. This was described in details in Ref. [50, 53]. 

Scope of the present script is the regeneration for the lattice QCD calculations [47-49]. In the present 
work, we tackle this problem through comparison with various effective models. In doing this, we have 
modified LSM and present systematic analysis for two order parameters. We have added to LSM the 
gluonic sector of the quasi-particle model. This is the essential original proposal of the present script. 
Thus, waiving details about PLSM itself is though as legitimated. But for a complete list of the PLSM 
parameters, the readers are kindly advised to consult [50, 53]. 
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B. Normalized net-strange condensate Aq^s{T) 



Fig. 2: The thermal dependence of Aq,s calculated from PLSM (solid curve), QLSM (dotted curve), PNJL (double- 
dotted curve) and HRG (dash-double-dotted curve) and compared with the lattice QCD calculations (solid circles) [47] 
and (open circles) [49] at mj/ms = 0.037. 


Another dimensionless quantity shows the difference between non-strange and strange condensates. 


A,,.(r) 


(qq) - ^{SS) 
m. 


{qq}o - 

m. 


0 


(56) 


where (qq) {{ss)) are non-strange (strange) condensates, and niq (rus) are non-strange (strange) masses. 
Using Ward identities and Gell-Mann-Oakes-Renner relation, expression (56) might be given in terms 
of pion and kaon masses and their decay constants [82]. Accordingly, the final results might be scaled 
but their thermal behavior remains unchanged. The lattice QCD calculations for Ag^s{T) (solid circles) 
[47] and (open circles) [49] are compared with the calculations from PLSM (solid curve), QLSM (dotted 
curve), PNJL (double-dotted curve) and HRG (dash-double-dotted curve) in Fig. 2. 

It is obvious that PLSM agrees with the lattice results [49] at low and also at high temperatures. Its 
ability to reproduce the other set of lattice results [48] is limited to the high temperatures. This might 
be originated in the difference between the two sets of lattice QCD simulations, section III A 1. The HRG 
model agrees well with this lattice calculations [48]. It is apparent that such agreement is limited to 
temperatures below the critical value due to the limited applicability of the HRG model. The remaining 
two models PNJL and QLSM show qualitative thermal behavior, as that from the other effective models 
and lattice calculations, which can be described by large plateau at low temperatures, around the critical 
temperature the values of Ag^s{T) decrease rapidly and at high temperature, Ag^s(T) vanishes but very 
slowly. Both models are closer to [49] rather than to [48] 

Both Figs. 1 and 2 show the PNJL model and HRG model describe much better the LQCD data for the 
magnetization and normalized net strange condensate, respectively, than for the chiral condensate. One 
should bear in mind that the magnetizations have been simulated in a different lattice that the one for the 
net strange condensate. Unfortunately, both quantities are not available from the same lattice simulation. 
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C. QCD chiral phase-diagram 


For mapping out the QCD chiral phase-diagram, various approaches are available. From PLSM and PNJL, 
as they possess two order parameters; one for strange and one for non-strange chiral condensates, hints 
about QCD chiral phase transition can be analysed. Furthermore, PLSM and PNJL possess deconlinement 
order parameter because of the Polyakov loop potential. Therefore, from strange and non-strange chiral 
condensates, a dimensionless quantity reflecting the difference between both condensates, Ag_s(T), can 
be deduced as function of temperature at fixed baryon chemical potential. Apparently, this signals the 
QCD chiral phase-transition. At the same value of baryon chemical potential, we can also deduce the 
deconfinement order-parameter as function of temperature. At a fixed baryon chemical potential, the 
thermal dependence of these two quantities intersect with each other at a characterizing point representing 
the phase transition. When repeating this procedure at different values of the baryon chemical potentials, 
we get a set of points representing the QCD phase-diagram. The results are given in Fig. 3, as solid curve 
for PLSM and dotted curve for PNJL. 

For the QCD chiral phase-diagram from QLSM, we implement another method. As no Polyakov loop 
potential is included, the QCD chiral phase-diagram is characterized by the higher-order moments of par¬ 
ticle multiplicity [31], which are assumed to highlight various types of fluctuations in T and Therefore, 
we utilize the possible fluctuations accompanying normalized second-order moment [31] in mapping out 
the QCD chiral phase-transition. The problematic of determining pseudocritical temperature from the 
second moment has been discussed in Ref. [83]. Accordingly, we observe that the peaks corresponding 
to different temperatures are conjectured to be characterized by different values of the baryon chemical 
potentials, where the QCD chiral phase-transition is conjectured to occur. We analyse this dependence 
at different values of the temperature T. Then, we follow the scheme to determine T and which is 
characterized by maximum m 2 /^^, where m 2 is the second-order moment of the particle multiplicity. The 
results are illustrated in Fig. 3, as dash-dotted curve. 

For the HRG model, we map out the QCD chiral phase-diagram by utilizing the quark-antiquark con¬ 
densate as order parameter [45]. It is assumed that the thermal dependence of the quark-antiquark 
condensate remains finite in the hadronic phase, but vanishes at temperature higher than the critical 
chiral-temperature. The results are given in Fig. 3, as well, as double-dotted curve. 

We can now shortly summarize the methods implemented to determine the pseuodcritical temperatures: 

— PLSM and PNJL: due to chiral and deconfinement phase transitions for light and strange quarks, 
Ag^s(T) is determined as function of temperature at a fixed baryon chemical potential, jjL. This signals 
the QCD chiral phase-transition. At the same value of /r, the deconfinement order-parameter can be 
studied as function of temperature, as well. Then, the thermal dependence of these two quantities 
is conjectured to intersect with each other at a characterizing point. When repeating this procedure 
for different values of /i, a set of points of pseudocritical temperatures and fi can be deduced. 

— QLSM: the normalized second-order moment of particle multiplicity is implemented in mapping out 
the QCD chiral phase-transition. Peaks corresponding to different temperatures are conjectured to 
be characterized by different values of T-^ and fj,. 

— HRG: the quark-antiquark condensates are implemented as order parameters. At vanishing and finite 
/i, the thermal dependence of the quark-antiquark condensate remains finite in the hadronic phase 
and vanishes at temperature higher than the critical chiral-temperature, T^. 

We observe that the chiral boundary from PLSM (solid curve) is positioned within the upper band of 
the lattice QCD calculations [84, 85] and agrees well with the freeze-out results deduced from the STAR 
BES measurements (symbols) [86]. The temperatures calculated from the HRG model by using the 
quark-antiquark condensate as the order parameter (double-dotted curve) [45] is higher than the chiral 
temperatures from the PLSM and the freeze-out temperatures calculated in the lattice QCD (band) and 
from the STAR BES measurements (symbols). Despite this difference, the corresponding T-fi sets are 
very similar to that of the PLSM. The results from PNJL and QLSM are higher than that from the HRG 
model. 
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Fig. 3: The PLSM T-^ chiral phase-diagram (lines with points), with which the freeze-out parameters deduced from 
lattice the QCD calculations [84, 85] (band) and that from different the thermal models [86, 87] (symbols) are compared. 


IV. CONCLUSIONS AND OUTLOOK 


In the present work, we report on a systematic comparison between PLSM, PNJL, QLSM and HRG 
in generating the chiral quantities, order parameter M (T) and normalized net strange and non-strange 
condensates Ag_s(T). Furthermore, we confront the results deduced from the four effective models to 
the recent lattice QCD calculations in order to distinguish between the models and to interpret the first- 
principle lattice QCD calculations. 

For the order parameter M(T), the best agreement is found with PLSM, while the recent lattice QCD 
[48] lay on top of all curves. This might be understood from the lattice configurations and the actions 
implemented in the simulations. The other effective models lay below the two sets of the lattice calcu¬ 
lations. We notice that the effective PLSM, QLSM, PNJL and HRC and the different sets of the lattice 
calculations have different critical temperatures. Tab. 11. 

For the normalized net strange and non-strange condensates Aq^s(T), PLSM again gives an excellent 
agreement with the lattice results [49] at low and high temperatures. But its ability to reproduce the 
lattice simulations [48] is limited to high temperatures. This might be originated in the difference between 
the two sets of lattice QCD simulations. Furthermore, we find that the HRC model agrees well with 
this the lattice QCD calculations [48]. It is apparent that this is restricted to temperatures below the 
critical value. The effective models PNJL and QLSM show the same qualitative thermal behavior. There 
is a large plateau at low temperatures. Around the critical temperature the values of Ag^siT) decrease, 
rapidly. At high temperature, Ag^siT) vanishes but very slowly. The effective models PNJL and QLSM 
are closer to [49] rather than to [48]. 

In light of this, we conclude that the PLSM reproduces M(T) and Ag^s{T), well. The HRC model is able 
to reproduce Ag^siT), while PNJL and QLSM seem to fail. These features and differences are present in 
the chiral phase-diagram. Fig. 3, as well. 

In section HIC, we have introduced the various order parameters used in the different models in order 
to deduce T and /r of the QCD chiral phase-transition. The strange and non-strange chiral condensates 
and the Polyakov loop potentials are utilized in PLSM and PNJL. The thermal dependence of these 
two quantities are assumed to intersect with each other at a characterizing point representing the QCD 
chiral phase-transition. For QLSM, no Polyakov loop potential is included in, therefore, the chiral phase- 
diagram is characterized by the higher-order moments of the particle multiplicity. The possible fluctuations 
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accompanying the normalized second-order moment are assumed to map out the QCD chiral phase- 
transition. For the HRG model, we utilize the quark-antiquark condensates as order parameter. 

Again, we find that the PLSM chiral boundary (solid curve) is located within the upper band of the 
lattice QCD calculations and agrees well with the freeze-out results deduced from the experiments and 
the thermal models (symbols). It is obvious that the chiral temperature calculated from the HRG model 
is larger than that from the PLSM. This is also larger than the freeze-out temperatures calculated in the 
lattice QCD (band) and from the experiments and the thermal models (symbols). Despite this difference, 
the corresponding T and fjL sets are very similar to that from the PLSM. This might be explained as follows. 
The T and ^ are calculated using different order parameters; in HRG vanishing quark-antiquark condensate 
but in the PLSM crossing (equalling) the chiral condensates and the Polyakov loop potential. The latter 
assumes that the two phase transitions; the chiral and the deconfinement, occur at the same temperature. 
The earlier deals with the chiral phase-transition independent on the confinement-deconfinement one. 

The results from the two model PNJL and QLSM show the same qualitative behavior. The chiral tem¬ 
peratures are higher than that from the PLSM and HRG. This might be interpreted due the heavy quark 
masses implemented in both models. 

Any model comparison with lattice results should span as much as possible of the parameter space. Even 
with the narrow parameter space explored in the present paper, we would like to highlight that the results 
are limited this. But, with reference to previous work [50], the parameters alone are not able to explain 
the diversity with the results in this study. We have to attack essential components of LSM and integrate 
gluonic sector taken from the quasi-particle models. 
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